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Abstract 

In  this  paper  plane  elastic  curves  are  revisited  from  a  viewpoint 
that  emphasizes  curvature  properties  of  these  curves.  The  family  of 
elastic  curves  is  considered  in  dependence  of  a  tension  parameter  a  and 
the  squared  global  curvature  maximum  k^.  It  is  shown  that  for  any 
elastic  curve  k^  is  bigger  than  the  tension  parameter  a.  A  curvature 
analysis  of  the  fundamental  forms  of  the  elastic  curves  is  presented. 
A  formula  is  established  that  gives  the  maximum  turning  angle  of  an 
elastica  as  a  function  depending  on  k^  and  a.  Finally,  it  is  shown 
that  an  elastic  curve  can  be  represented  as  a  linear  combination  of 
its  curvature,  arc  length  and  energy  function  and  that  any  curve  with 
this  property  is  an  elastica. 

1      Introduction 

The  search  for  Smoothing  Algorithms  in  CAGD  has  induced  research  on 
curves  and  surfaces  which  minimize  functionals  with  geometrical  or  physical 
meaning.  The  classical  example  are  elastica  which  describe  the  shape  of 
elastic  materials.  Plane  elastica  can  be  defined  as  the  extremals  (critical 
points)  of  the  variational  problem 

tL 

I    k2(s)  +  cr  ds  — +      min 
Jo 

where  k  denotes  the  curvature,  5  the  arc  length  of  a  plane  curve  x  and  cr  is 
constant.  The  total  length  L  of  i  is  considered  to  be  variable.  The  set  of 
comparison  curves  of  the  problem  is  the  set  of  all  C°°  curves  in  the  plane 
with  fixed  endpoints  and  fixed  tangent  directions  at  these  points. 


Lee  and  Forsythe  found  (see  [6])  that  the  curvature  function  of  an  elastica 
satisfies  the  differential  equation 

k"{s)  +  ^k3(s)  +  Kk{s)  =  0,     KeR.  (1) 

We  begin  this  paper  with  a  short  derivation  of  this  result  and  a  proof  that  this 
equation  is  also  sufficient  for  a  curve  to  be  an  elastica.  We  will  then  show  that 
for  any  elastica  the  square  of  the  curvature  function  has  a  global  maximum 
k^  which  is  bigger  than  the  tension  parameter  a.  The  fundamental  forms  of 
elastic  curves  first  described  by  Euler  are  characterized  by  the  values  of  the 
parameters  /c^  and  a.  A  curvature  analysis  of  these  curves  is  presented. 

In  section  3  various  formulas  expressing  properties  of  elastica  in  terms  of 
k^  and  a  are  given.  This  includes  a  formula  that  gives  the  maximum  turning 
angle  of  an  elastica  as  a  function  depending  on  K2m  and  a.  The  explicit 
formula  for  the  maximum  turning  angle  allows  to  determine  an  upper  bound 
for  n2m  in  the  case  that  an  inflectional  elastica  has  a  turning  angle  bigger 
than  7r.  As  the  main  result  of  this  section  we  show  that  a  plane  elastica  has  a 
representation  as  a  linear  combination  of  its  curvature,  arc  length  and  energy 
function  E(s)  =  f03  k2(s)  ds.  This  property  characterizes  the  class  of  plane 
elastica  and  can  be  used  for  speeding  up  the  computation  of  interpolating 
elastica. 


2      The  curvature  of  plane  elastica 

Given  two  points  P,  Q  €  R2  and  two  unit  vectors  V  €  TpR2,W  G  TqR2 
then  M  denotes  the  set 

M  :={x:  [0,L]-+R2  :        L  <E  R+,i  €  C°°[0,  L],  \x'(s)\  =  1  for  s  €  [0,1], 

*(0)  =  P,  x{L)  =  Q,  *'(0)  =  V,  x\L)  =  W). 

We  consider  the  problem  of  minimizing  the  functional 

tL 


E  +  vL:=  I    K2(s)  +  <rds 
Jo 


among  all  curves  of  M  where  k  denotes  the  curvature  and  5  the  arc  length 
parameter  of  a  plane  curve  x.  The  integral  E  =  /0  k2{s)  ds  is  called  the  en- 
ergy of  x  on  [0,  L).  a  denotes  a  constant  which  we  call  the  tension  parameter. 
The  total  length  L  is  considered  to  be  variable.  As  the  functional  considered 
is  translation  and  rotation  invariant  we  assume  in  this  section  that  P  =  (0,0) 
and  V  =  (1,0). 


For  x  G  M  the  tangent  vector  T  of  x  is  given  by 

T(s)  =  (cos(tf(s)),  sin(*(s))) 

where  the  function  ^  with  ty(s)  :=  f0s  k(s)  ds  gives  the  turning  angle  of  x. 
Using  this  notation  the  variational  problem  can  be  written  as 


rL 

min       /    W^s)  +  a  ds 

*€C~[0,L]./o 


under  the  constraint 


Jo    \s\n*(s)J  CS      V' 


<(s) 
where  the  admissible  functions  ty  G  C°°[0,  L]  are  subject  to  the  constraints: 

#(0)  =  0,      $(L)  =  ip     with     ( cos  V>,  sin 0 )  =  W'T. 

According  to  Lagrange's  multiplier  rule  the  differential  equation  of  this 
problem  is  the  Euler  equation  for  the  integrand 

F{y,V,s)  =  #'2(s)  +  (T  +  Acosty(s)  +  /isinty(s), 

for  some  constants  A,/i  £  R,  i.e. 

\K(5)  =  -^  sin  V{s)  +  I  cos  *(s).  (2) 

(see  e.g.  in  [3], [6]).  By  introducing  the  constants 

A  =  2a  cos  4>     fi  =  2a  sin  <f> 
we  rewrite  the  Euler  equation  in  the  form 

V"  =  k'  =  -asm{V-<f>).  (3) 

Multiplying  (3)  by  2^'  and  integrating  yields 

/c2  =  2acos(*  -<}>)  + A  (4) 

where  A  €  R  denotes  an  integration  constant.   Equation  (4)  has  been  used 
to  define  elastic  curves  in  the  plane  in  the  classic  literature  (see  [7]). 

In  order  to  determine  the  constant  A  in  terms  of  the  tension  parameter 
<j,  we  consider  the  boundary  condition 

F(¥(L),¥'(L),L)  -  -^F(*  (s),  *'(*),  s)\l*'(L)  =  0 


that  must  be  satisfied  by  the  extremal.  This  condition  is  implied  by  the  fact 
that  the  total  length  L  of  the  curve  is  variable  in  the  variation  (see  e.g.  [2], 
p.  571).  Thus, 

*"(/,)  =  k2(L)  =  A  cos  Vf(L)  +  ft  sin  *(I)  +  a.  (5) 

Comparing  (5)  with  (4)  shows  that 

A  =  a. 

Therefore  we  give  the  following  definition  of  an  elastic  curve  under  tension. 

Definition  1  An  arc  length  parametrized  plane  curve  x  with  curvature  func- 
tion k  is  called  elastica  (or  elastic  curve)  with  tension  parameter  a,  if  for 
some  a,(^ER 

k2  =  2acos(*  -  4>)  +  o  (6) 

where  ty  denotes  the  function  ty(s)  :=  /0*  K(s)ds. 

Lee  and  Forsythe  showed  in  [6]  that  the  Euler  equation  (2)  implies  the 
differential  equation  (1)  for  the  curvature  function  /c.  We  will  now  show  that 
(1)  is  in  fact  equivalent  to  (6). 

Theorem  2  If  k  €  C2(R)  and  $(s)  :=  /0S  K{s)ds,  then 

k2  =  2acos($  -  <f>)  +  0- 

holds  for  some  constants  a,  ^  and  a  if  and  only  if 

k"  +  \k3  -  \vk  =  0.  (7) 

Proof:  (i)  Differentiating  (6)  yields  the  Euler  equation  in  the  form  of 
(3).  Differentiating  the  Euler  equation  gives 

k"  =  —  acos(ty  —  4>)k. 

Substituting  the  term  cos(ty  —  <p)  in  the  above  equation  according  to  (6)  we 
obtain  (7). 

(ii)  Differential  equation  (7)  can  be  integrated  to  the  first  oder  differential 
equation 

(K')2  =  C-(\/4)a2-(\/4)(K2-a)2.  (8) 

Note,  that  for  a  real  solution  of  (8)  it  is  necessary  that 

a2  :=C-(l/4)<72>0 


and 

("2  ~  °Y 
4a2 
The  function  6  defined  by 


<  1. 


K     —  <J 


0  :=(-l)n  arccos(-5?-r) 

with  n  =  1  for  «'  >  0  and  n  =  2  for  /c'  <  0  obviously  obeys  the  relation 

k2  =  2\a\cos0  +  a.  (9) 

It  remains  to  show  that  0  is  an  integral  of  k. 

(8)  together  with  the  definition  of  a2  and  (9)  implies 

(*')2  =  a2sin20. 

Since  sin  0  is  positive  (negative)  if  k'  is  negative  (positive),  we  obtain 

K'=-\a\sm6.  (10) 

Differentiating  (9)  and  substituting  k'  according  to  (10)  yields 

0'  =  K. 

Since  ^  and  6  are  both  integrals  of  k  there  is  a  constant  4>  such  that 

<a  =  e  +  <t>.  For  vp 

k2  =  2|a|cos(^  -  <£)  +  <7 
holds  because  of  (9). 


The  squared  curvature  function  /c2  of  an  elastic  curve  has  a  global  maxi- 
mum, even  if  /c2  is  extended  to  the  whole  real  line. 

Lemma  3  If  k  is  a  solution  of  (7)  on  R,  then  k2  has  a  global  maximum. 

Proof:  Since  the  function  (1/2)(k3  —  an)  is  of  the  class  C](R)  the  solu- 
tions of  (7)  can  be  extended  to  the  whole  real  line. 

To  show  that  k2  has  a  global  maximum  we  observe  first  that  (3)  implies 
the  existence  of  a  local  extremum  of  k.  This  is  because  the  assumption 
k'(s)  =  ^"(s)  /  0  for  all  s  £  R  means  that  ty  is  convex  or  concave  and 
therefore  unbounded  while  according  to  (6)  k'  has  no  zeros  only  if  ty  is 
bounded.  Furthermore  from  (3)  and  (4)  it  is  obvious  that  any  local  extremum 
of  k  is  a  global  extremum  of  k2. 


We  assume  now  that  k2  has  no  global  maximum  on  R.  In  this  situation 
k'  has  exactly  one  zero  smtn  on  R  and  k2  takes  its  global  minimum  in  this 
point.  Note  that  k  is  non-zero  for  any  s  ^  smtn  because  a  zero  of  k  at  a 
point  s  ^  smtn  would  imply  that  /c(smtn)  =  0  and  k'(s)  =  0  for  some  point  s 
between  s  and  smtn.  Therefore  (k2)'  has  no  zero  besides  smtn  and  k2  is  strictly 
monotone  increasing  on  the  right  of  smtn.  Hence  k  is  monotone  increasing 
resp.  decreasing  on  the  right  of  smin  if  k  has  positive  resp.  negative  values 
on  the  right  of  smfn.  The  formula 


*($)=    f       K{s)ds  +  V{smin) 
«'*min 


yields  that  ty  is  in  any  case  unbounded.  This  is  a  contradiction  to  the  as- 
sumption because  (4)  implies  the  existence  of  global  maxima  for  k2  if  $  is 
unbounded. 

D 

We  now  give  the  relation  between  the  global  maximum  k2^  of  k2  and  the 
tension  parameter  a  and  express  k  in  terms  of  elliptic  functions. 

Theorem  4  Let  k  €  C2(R)  be  a  solution  of  the  differential  equation  (7)  with 
a  global  maximum  k2^  =fi  0  of  k2  on  R.    Then  the  following  statements  hold: 

(i)  «l  >  ° 

(ii)  k  has  a  zero  if  and  only  if  k^  >  2a.  In  this  case  k  is  given  by 


k(s)  =  Km  cn(y/(K2m-a)/2(s  -  sm)  \  k2)  (11) 

with  the  parameter 


k2  = 


(Hi)   2a  —  k2^  <  K2  <  k2^  for  k2^  <  2a.  In  this  case  k  is  given  by 

k{s)  =  Kmdn(Km{s  -  sm)/2  |  — ).  (12) 

Proof:  As  a  global  maximum  point  sm  of  k2  is  also  a  zero  of  k'  (8)  implies 

C=-K4m-  -CK2 

4  m      2      m 
Hence  (8)  takes  the  form 

(*')'  =  j(«m"  «2)(«2  +  4  -2^).  (13) 


As  all  quantities  in  (13)  are  real  (i)  follows  from  the  fact  that  the  term 
(Km  ~  k2)  *s  always  non-negative  and  therefore  the  term  (k2  +  K2m  —  2a)  has 
also  to  be  non-negative  for  any  k2. 

Using  the  same  argument  we  observe  that  k2(s)  =  0  for  some  5  implies 
that  k2^  >  2<r  while  for  k^  <  2a   k2  has  to  be  greater  or  equal  2<r  —  «2n. 

To  express  k  in  terms  of  elliptic  functions  we  procede  as  follows.  In  the 
case  that  /c2^  >  2a  we  substitute  z2  —  (/c^  —  «2)/«2n  in  (13)  and  obtain  for  z 
the  differential  equation 


M^Jod-.rXl-z'Xl-fcV)  (14) 


where 


k2  = 


2«-^)- 

(14)  is  the  differential  equation  of  Jacobi's  function  sn=sn(u)  for  the  argu- 
ment u  =  -J^i^  —  a)s  (see  [l],p.  114).  We  therefore  obtain 

z*(s)  =  sn2(y/(Kll-c)/2(s-3n)\k2). 
The  relation  k2  =  n2m{\  —  z2)  then  gives 

k2(s)  =  K2mcn2(J(K2m-a)/2(s  -  sn)  \  k2)  (15) 

(see  [8],p.l6).  Since  k  is  differentiable,  (15)  implies  (11). 

In  the  case  that  K2m  <  2a  we  substitute  z2  =  (/c^  —  k2)/2(k2ti  —  a)  in  (14) 
and  obtain  for  z  the  differential  equation 


(*T  =  1,4(1  -  z')(l  -  IV) 


where 

Hence 

The  relation 

yields 


J_      2(K2m  -  a) 
k2  '  -2 


*m 


z2(5)  =  sn2(Km(s-sm)/2|/2). 

K2  =  K2m(l-l2Z2) 


K2(s)  =  K2mdn2(Km(s-sm)/2\l2) 
(see  [8],p.l6)  which  implies  (13)  again  because  of  the  differentiability  of  n. 


Theorem  2  implies  that  the  curvature  function  of  an  elastic  curve  extends 
to  a  periodic  function  on  R.  In  the  case  that  «£,  >  2a  this  periodic  function  is 
symmetric  with  respect  to  any  zero  of  its  derivative  and  antisymmetric  with 
respect  to  any  zero.  According  to  Love  [7]  this  situation  is  called  inflectional 
because  the  extension  of  x  has  turning  points.  In  the  case  that  K2m  <  2a  the 
extension  of  k  has  no  zeros  but  is  still  symmetric  with  respect  to  any  zero 
of  its  derivative.  The  local  extrema  in  this  case  are  /cm  and  ±(2cr  —  k^). 
The  situation  is  illustrated  in  figure  1  -  figure  8  where  for  /cm  =  1  various 
curvature  functions  and  the  corresponding  elastica  are  shown.  If  only  non- 
negative  tension  values  are  considered  a  is  an  element  of  [0,1].  With  a 
increasing  from  0  to  1  the  curvature  function  changes  continuously  from  a 
lemniscate  function  (a  =  0)  to  a  constant  (a  =  1).  A  classification  of  the 
different  forms  of  elastica  has  been  given  by  Euler  (see  [5],  [7]). 


Figure  1.  a  =  0;  /c(s)  =  (^  coslemn(/cm(5  -  sm)/2). 


1C  12 


Figure  2.  a  =  0.3;  *(*)  =  Km  cnfv/^  -  ct)/2(*  -  sm)  \  k2) 


•M  -10  $ 


Figure  3.  a  —  0.35;  k(s)  as  in  figure  2. 


Figure  4.  a  =  0.394757217;  «(*)  as  in  figure  3. 
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Figure  3.  a  =  0.49;  #c(«)  =  *m  cn^/c^  -  a)/2(*  -  «m)  |  k7). 


-15  -10  -5  0  J  10  I! 


Figure  4.  a  =  0.5;  *(s)  =  Km  sech(/cm(5  -  sm)/2). 


•15  -1C  -5  C  3  10  15 


Figure  5.  c  =  0.62;  «(s)  =  ^  dn(*m(5  -  5m)/2  |  £). 


-15  -10  -J  0  J  JO  U 


Figure  6.  a  =  1;  k(s)  =  k^. 
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3      A  characterizing  representation  theorem 

First  we  summarize  the  main  formulas  for  elastic  curves.  In  order  to  avoid 
repetitions  we  use  the  following  convention:  if  a  curve  x  :  [0,  L]  — ►  R2  is  an 
elastica  ,  then  k  denotes  the  (analytic)  curvature  function  of  x,  k  denotes  the 
extension  of  k  on  R,  k^  the  global  maximum  of  k2  and  sm  is  a  number  with 
k2(sm)  =  k2^.  Furthermore  ^  :  [0,  L)  — ♦  R  is  the  function  ty(s)  :=  /0*  k(s)  ds, 
f"  denotes  the  extension  of  ty  on  R  and  <p  :=  $(sm).  The  energy  E  of  x  is 
the  function  E(s)  =  /0*  k2(s)  ds.  For  convenience  we  also  use  the  notation 
C(rp)  =  (cosV>,sinV0.  i?  is  the  angle  in  [0,2tt[  with  x'(0)  =  C{ti). 

Theorem  5  For  a  plane  elastic  curve  x  with  tension  parameter  a  the  fol- 
lowing relations  hold: 

(i)   k'(s)  =  -}(!&-  a)  sin(*  -tp), 

(ii)    K2(s)  =  (K2m-o)cOs(*(s)-v)  +  Cr, 

(iii)  k(s)  -  «(0)  =  -|(i£  -  a)  <  C'(<p  +  r?), x(s)  -  x(0)  >. 
(iv)  E(s)  =  (4  -  a)  <  C(y>  +  *), *(a)  -  x(0)  >  -ras, 

Proof:  Without  loss  of  generality  we  can  assume  that  a  in  (4)  is  positive. 
From  lemma  3  we  know  that  k2  has  a  global  maximum  k^.  Since  for  positive 
a  the  maximum  occurs  iff  cos(ty(sm)  —  <j>)  =  1,  it  follows  that 

2a  =  K2m-a 

and 

<£  =  tf(sm)  +  2*?r. 

Therefore  (3)  and  (4)  imply  (i)  and  (ii).  To  obtain  (iii)  we  rewrite  (i)  as 

K'(s)  =  --{n2m  -  (r)(sin(tf  +  t?)  cos(<^  +  <&)  -  cos(tf  +  d)  sin(v?  +  d)) 

and  integrate  using 

x'(5)  =  (cos(*(5)),sin(*(s))). 

To  verify  (iv)  we  write  (ii)  in  the  form 

K2(s)  =  (n2m  -  <r)(cos(tf  +  ti)  cos(y?  +  t?)  +  sin(*  +  d)  sm(ip  +  #))  +  <r 
and  integrate. 
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Since  the  elliptic  function  en  in  Theorem  3  has  the  symmetry  properties 
of  a  sine  wave  while  dn  is  positive,  the  turning  angle  of  an  elastic  curve  is 
bounded  in  the  inflectional  case  a  <  ^ac^  and  unbounded  otherwise.  Formula 
(ii)  of  Theorem  3  can  be  used  to  determine  the  maximum  turning  angle  of 
an  inflectional  elastic  curve. 

Corollary  6   The  maximum  turning  angle  ^max  :=  maxj6[(i,+^  |^(s)|  of  an 
elastica  x  with  period  T  and  tension  parameter  a  <  \K2Ji  is  given  by 


a 


*maX  =  2arccos( -).  (16) 


a  —  k2 


Ifa>  \K2m  /  0  then 

|*(<  +  :r)-*(t)|  =  27r 

for  any  /£R. 

Proof:  Let  s0  be  a  zero  of  k.  Then  it  follows  from  the  symmetry  prop- 
erties of  k  in  the  inflectional  case  that 

rao+T/4 

*mar  =  2|    /  R(s)  ds\  =  2\V  (s0  +  T/4)  -  ¥(«l)|. 

J  so 

With  (ii)  we  get 

Vmax  =  2  arccos  — )  -\-ip-  arccos(  — )  -  tp\ 

and  therefore  (16). 

In  the  non-inflectional  case  k  is  a  positive  periodic  with  the  period  T  = 
4A7|Km|  where  K  denotes  the  complete  elliptic  integral  of  the  first  kind.  For 
t  =  sm  one  obtains 

|*(sm  +  T)  -  tf  (sro)|  =  4  /    dn(u|/2)  du  =  4arcsin(sn(/0)  =  2*. 

Jo 


For  an  inflectional  elastica  x  formula  (16)  implies  an  upper  bound  for  k^ 
if  the  absolute  value  of  the  oriented  angle  between  x'(0)  and  x'(L)  is  bigger 
than  7T. 

Corollary  7  For  an  inflectional  elastica  x  with  tension  parameter  a  and 

2<*k<l ! . 

a  cos(rp/2) 
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Proof:  It  follows  from  (16)  that  a  turning  angle  bigger  than  tt  can  only 
happen  for  positive  a.  The  inflectional  nature  of  x  implies  therefore  the  left 
inequality.  The  right  inequality  follows  from  the  fact  that  %l>  has  to  be  less 
or  equal  <Pmar  which  is  given  by  (16). 

D 

As  the  main  result  of  this  paragraph  we  show  that  an  elastica  can  be 
represented  as  a  linear  combination  of  its  curvature,  energy  and  arc  length. 

Theorem  8  If  x  is  an  elastica  with  tension  parameter  a  ^  k^  then 

{)~  Kl-*\-cos(<p  +  d)     sin(y>  +  t?)A     E(s)  -  as     J+'WVn 
where  cos  \p  —  (kI  —  a)/(K2n  —  a)  and  sin  \p  —  2k'0I{k2ti  —  a). 
Proof:  A  plane  curve  x  £  C°°[0,  L]  has  a  representation 

x(s)-x(0)    =     ['  C(V{s)-  if  +  <p  +  tf)  ds 
Jo 

=    C{<p  +  d)  I  cos(*(s)  -  ip)  ds 

Jo 

+    C'Up  +  ti)  [*  sm( V(s)-<p)ds  . 
Jo 

Applying  (iii)  and  (iv)  of  Theorem  3  yields 
x(s)  -  x(0)  =  -jl—  (C(if  +  d)(E(s)  -  as)  -  C'{tp  -r  *)2(k(s)  -  /c(0))) 


Km"* 


which  is  equivalent  to  (17). 

The  formulas  for  (p  follow  from  (i)  and  (ii)  for  6  =  0. 


Formula  (17)  provides  an  explicit  representation  of  an  elastica  in  terms 
of  its  curvature  function  if  k2^  ^  a.  Note  that  in  the  case  /c^  =  a  the  elastic 
curve  is  a  circle  of  radius  l/|«m|  &s  shown  in  figure  8. 

Finally  we  show  that  elastica  are  essentially  the  only  curves  in  the  plane 
that  have  a  representation  of  the  form  (17). 

Theorem  9  Let  k  be  an  arbitrary  C2(R)  function  with  a  global  maximum 
k2^  of  k2 ,  a  a  real  number  smaller  than  n2m  and  E(s)  :=  /0'  n2(s)ds.  An  arc 
length  parametrized  curve  x  where  x(s)  is  given  by  (17)  is  an  elastica  with 
curvature  function  k  and  tension  parameter  a. 
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Proof:  We  assume  the  case  <p  -f  i?  =  tt/2  which  can  always  be  achieved 
by  applying  a  rotation  to  x.  Then 

x  =  — (2k',  k2  -  a),     x"  =  — - (2k",  2kk'). 

x  is  arc  length  parametrized  if  and  only  if  x'  =  1,  i.e. 

(2k')2  =  (£  -  of  -  (k2  -  a)2.  (18) 

Differentiating  (18)  yields 

2k k"  +  k'k(k2  —  a)  =  0. 

Therefore  either  k  satisfies  (7)  or  k  is  constant.  For  an  arc  length 
parametrized  curve  the  determinant  [x1,  x")  is  curvature.  Here  we  have 

If  k  is  constant  the  curvature  of  x  is  zero,  hence  x  is  a  trivial  elastic  curve. 
If  k  satisfies  (7)  we  substitute  in  (19)  k"  according  to  (7)  and  k'2  according 
to  (18).  This  yields 

[X,X       =  K. 


The  representation  formula  (17)  of  an  elastic  curve  is  extremely  useful  for 
the  computation  of  interpolating  elastica.  One  reason  for  this  is  that  (17) 
involves  no  trigonometric  functions  and  fewer  integrations  than  the  standard 
representation  based  on  the  formula  x'  =  (cos  ty,  sin  ^).  Therefore  (17)  is 
less  expensive  to  evaluate.  Furthermore  (17)  can  be  used  to  find  piecewise 
polynomial  approximations  of  elastic  curves  based  on  a  spline  approximation 
of  the  curvature  function.  The  author  has  established  polynomial  splines 
which  approximate  the  curvature  functions  of  plane  elastica  with  high  ac- 
curacy. Using  these  piecewise  polynomial  curvature  functions  together  with 
the  new  representation  (17)  one  obtains  polynomial  spline  approximations  of 
the  elastica  itself.  These  approximations  will  be  discussed  in  [4]. 
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